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Abstract. This study is the first known attempt to search for 
gamma-ray burst repeaters combining data from gamma-ray 
experiments flying on board different satellites and making use 
of information derived from the bursts detected simultaneously 
by all the experiments. The proposed method is suitable to cor- 
relate GRB data provided by experiments that overlap partially 
or totally in time. 

As an application of this method we have correlated the po- 
sitions of 57 gamma-ray bursts observed by WATCH/GRANAT 
and WATCH/EURECA with 1905 bursts detected by BATSE. 
Comparing the so-called "added correlation" between the 
WATCH and BATSE bursts with that obtained with simulated 
WATCH catalogues, we conclude that there is no indication of 
recurrent activity of WATCH bursts in the BATSE sample. We 
derive an upper limit of 15.8%, with a confidence level of 94%, 
for the number of WATCH gamma-ray bursts that could repre- 
sent a population of repeaters in the BATSE sample. 
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1. Introduction 

Despite the advances carried out so far, the origin of the 
gamma-ray bursts (hereafter GRBs) remains unknown. The 
identification of absorption lines in the optical spectrum of 
GRB 970508 strongly supports models arising from sources at 
cosmological distances (Metzger et al. 1997), but there is still 
a lack of knowledge on the mechanisms originating these enig- 
matic phenomena. One of the most important clues that could 
clarify the nature of the GRBs would be the detection of a re- 
peater behaviour. 

Initial studies showed an apparent evidence of repetition for 
the BATSE IB catalogue (Quashnock and Lamb 1993), sug- 
gesting that it would be possible to have an excess of pairs of 
GRBs clustered in both time and space (Wang and Lingenfel- 
ter 1995). This fact was not confirmed by the work carried out 
using the BATSE 2B catalogue (Brainerd et al. 1995), although 
other studies provided marginal evidence for both temporal and 
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angular clustering (Petrosian and Efron 1995). Analyses based 
on autocorrelations with data from the BATSE 3B catalogue 
did not find any evidence of repetition (Bennett and Rhie 1996) 
and have imposed several constraints to the number of repeaters 
(Tegmark et al. 1996). Finally, recent studies confirm the lack 
of repetition in the 4B catalogue and lead to an upper limit to 
the repetition rate of 0.04 burst source -1 yr _1 (Hakkila et al. 
1997). 

The BATSE 4B catalogue was obtained by the BATSE 
experiment on board the CGRO satellite and contains 1637 
GRBs detected from April 1991 to August 1996 (Paciesas et 
al. 1998). The BATSE experiment consists of eight identical 
detector modules, placed at the corners of the CGRO spacecraft 
and covering energy channels from ~ 25 keV to ~ 2 MeV. It 
provides error boxes with a minimum radius of 1.6° (lcr confi- 
dence level, Fishman et al. 1994). BATSE is detecting bursts at 
a rate of 0.8 bursts per day. The bursts are daily added to the so- 
called Current GRB Catalogue, which contains the BATSE 4B 
catalogue plus all bursts detected after August 1996. When this 
study was started, the catalogue contained 1905 sources; this 
sample constitutes the basis of the present work. 

The WATCH X-ray all-sky monitor is based on the rota- 
tion modulation principle (Lund 1986). The instrument has a 
circular field of view of 4 steradians and an effective area of ~ 
30 cm 2 (averaged over the field of view). Position sensitivity is 
achieved using the rotation collimator principle, with the colli- 
mator grids rotating with a frequency ui=l Hz. The phoswich 
detectors consist of interleaved scintillator-strips of Nal and 
Csl crystals. The geometric area of the scintillator is 95 cm 2 . 
Four units were mounted on board the Soviet GRANAT satel- 
lite in a tetrahedral configuration covering the whole sky, and 
one unit on board the European Space Agency EURECA space- 
craft. The total energy range is 8-80 keV, therefore overlap- 
ping with the lower BATSE energy band. WATCH/GRANAT 
detected bursts in 1990-94 and WATCH/EURECA in 1992-93, 
thus both experiments also overlapped in time with BATSE. 
One of the main advantages of WATCH was the capability of 
locating bursts with relatively small error boxes (3<r error radii 
with - 1°) (Brandt et al. 1990). WATCH/GRANAT detected 
47 GRBs in this period and WATCH/EURECA 12 (Castro- 
Tirado et al. 1994, Brandt et al. 1994, Sazonov et al. 1998). 
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Fig. 1. Error boxes for the 57 GRBs detected by WATCH, represented in galactic coordinates. The sample contains 45 GRBs 
detected by WATCH/GRANAT, 10 by WATCH/EURECA and two localized by both experiments at the same time. The typical 
radii of the error boxes are ~ 1°, with a 3a confidence level. 



Two GRBs (GRB 920814 and GRB 921022) were detected 
by both the WATCH/GRANAT and WATCH/EURECA exper- 
iments. Therefore, the sample of WATCH GRBs used in this 
study comprises 57 GRBs: 45 WATCH/GRANAT bursts, 10 
WATCH/EURECA bursts and the above-mentioned two GRBs. 
BATSE also detected 27 of them. Fig. [j] shows the sample of 
57 WATCH GRBs used in this study. 

The distribution of time amplitudes for GRBs shows two 
classes of bursts: a) durations shorter than ~2s and b) longer 
than ~ 2 s (Kouveliotou et al. 1993). It was noticed that the 
energy spectra of the short bursts were generally harder than 
those of the long ones (Kouveliotou et al. 1993, Lestrade et al. 
1993). 

The fraction of short events in the WATCH sample is 
smaller than that in the 4B catalogue. This fact can be justi- 
fied by at least three selection effects: 

i) The availability of WATCH for localizing sources is gov- 
erned by the rotation speed of the collimator grids (1 Hz). So, a 
source needs to be bright enough for at least one rotation of the 
modulation collimator in order to be localized, implying a burst 
duration longer than 1 s. In contrast, the BATSE experiment is 
able to detect bursts with durations as short as 64 ms. 

ii) The low energy band of the WATCH experiment (~8- 
20 KeV) is sensitive to the soft GRBs, below the BATSE lower 



limit (~25 KeV), which generally belong to the class of bursts 
with durations longer than 2 s. 

iii) On the other hand, since WATCH is about an or- 
der of magnitude less sensitive than the large-area detectors 
of BATSE, the WATCH catalogue contains bursts which are 
brighter than those in the BATSE sample. 

The above three reasons explain why the GRBs in the 
WATCH sample are longer, softer and brighter than the average 
BATSE 4B bursts. 

This study is the first known attempt to search for repeaters 
combining data of 7-ray experiments flying on board different 
satellites. The method proposed makes use of the so-called "si- 
multaneous bursts" and is suitable to correlate GRB data pro- 
vided by experiments that overlap partially or totally in time. 
In the future, this work could also be used to detect systematic 
pointing errors between different 7-ray experiments, allowing 
to improve their capability for locating GRBs. 

2. Method 

In this section we outline the methodology proposed to carry 
out the study. First, we exclude the simultaneous bursts (§ 2.1) 
and calculate the so-called "added correlation" function be- 
tween the WATCH and BATSE samples (§ 2.2). Afterwards, 
1500 WATCH catalogues are simulated (§ 2.3) in order to cal- 
culate the expected value of the "added correlation" function. 
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Then the distribution of the overlapping function for real and 
random overlaps is obtained (§ 2.4) and finally the probability 
of having different number of repeaters (§ 2.5) is found. 

We consider that there is a common source in both sam- 
ples when the emission of a repeater is detected at least twice, 
once by each experiment and the detections are separated in 
time. Thus, the same GRB detected simultaneously by both ex- 
periments is not considered as a common source. Our study is 
aimed at searching common sources detected by both WATCH 
and BATSE experiments. 

2.1. Simultaneous bursts 

The positions of 27 simultaneous bursts detected by WATCH 
and BATSE are in good agreement. If BATSE la error boxes 
are considered, there are 20 overlaps with WATCH 3<r boxes. 
Instead, if 3a error boxes are taken into account there is 
only one burst (GRB 920714) that does not overlap. These 
27 bursts were excluded from the BATSE sample of 1905 
sources, because they are obviously the same sources detected 
by WATCH. Therefore the sample was reduced to 1878 bursts. 
Nevertheless, the simultaneous bursts were considered in fur- 
ther calculations (§ 2.4), because they provide information 
on the overlapping expected for a repeater detected by both 
BATSE and WATCH. 

2.2. The "added correlation" estimate 

Recurrence, even in a single case, would be immediately obvi- 
ous if we had locations with no errors. However, the locations 
provided by BATSE and WATCH, while numerous, have inac- 
curacies and consequently a statistical analysis is required to 
demonstrate, or limit, the presence of common sources. If any 
repeater is present in both catalogues, an excess in the overlap 
between the error boxes of both catalogues would be expected. 
We define the overlapping function between the z-th WATCH 
and the j-th BATSE error boxes as the following integral over 
the galactic coordinates I and b: 




AFj J j Wi(l, k, b, bi)Bj(l, lj,b, bj) dfl if dy < <?, + <r, 



where A is a normalization factor computed in such a way 
that Cij remains between and 1. Fj is the BATSE exposure 
correction for the BATSE j-th burst. and bi,bj are the 
galactic coordinates of the centre of the i-th WATCH and the 
j-th BATSE error boxes, is the distance between the i- 
th WATCH and the j-th BATSE burst, and a u aj are the 3a 
radii of the i-th WATCH and the j-th BATSE error boxes. 
Bj(l, li, b, bi) is a Gaussian-like normalized probability distri- 
bution given by the following expression: 

Bj(l, lj,b, b,) = Hl ~ S) exp( {dj/ajf Ml - s) ) 

j 

with s = 0.9973, and dj the distance between the integra- 
tion point and the centre of the j-th BATSE burst; 



dj =arccos(sin(6j) sin(6) + cos(bj) cos(6) cos(Z — lj)) 

Wi(l,li,b,bi) is analogous to Bj(l,lj,b,bj) based on 
WATCH coordinates. Although we are aware that the errors of 
BATSE locations do not follow a single Gaussian distribution 
(see Briggs et al. 1998), we consider that, for our purposes, we 
can extend the Gaussian approximation from 1 a to 3a. This is 
a very appropriate and useful approximation which has been 
frequently used in the past (Fisher et al. 1987, Bennett and 
Rhie 1996), providing stringent upper limits on the 3B cata- 
logue (Tegmarket al. 1996). 

On the other hand, the error introduced in by consider- 
ing only overlaps between 3a error boxes, instead of assum- 
ing unlimited error boxes, is less than 0.1%, irrelevant for our 
final conclusions. In the approximation that a,-, « 60° and 
aj << 60° (which is quite accurate, since typical values are 
0-j ~ 1° and aj ~ a few degrees), dj approximately depends 

on dij like ~ exp '"^ f^Ji 3 , so it decreases rapidly when both 

probability distributions are not close to each other, dj pro- 
vides a measurement of whether both GRBs originated from 
the same source or not. Based on the former arguments, we 
define the "added correlation" C as follows: 

57 1878 

C^^dj 

3=1 i=l 

C is a parameter which is very sensitive to the presence 
of common sources in both catalogues. The larger the num- 
ber of common sources, the higher the value of C obtained. 
Our study is based on the comparison of the "added correla- 
tion" C calculated for the real WATCH catalogue (renamed as 
Cbw) with those obtained for 1500 WATCH simulated cata- 
logues (renamed as Cj,j — {1, 2, 1500}). C is the general- 
ization for two probability distributions (WATCH and BATSE) 
of the R statistics introduced by Tegmark et al. (1996). Cj is 
corrected by the BATSE and WATCH exposure maps, the first 
one is taken into account in the term Fj included in the defi- 
nition of c^, whereas the second one is considered to simulate 
the WATCH catalogues for which Cj are calculated. 

2.3. Simulation of WATCH catalogues 

Monte Carlo simulations of 1500 WATCH-like catalogues have 
been performed. They provided 1500 values for C called 
Cj,j = {1, 2, 1500}. In order to determine reliable val- 
ues for them, the exposure maps of the WATCH/GRANAT and 
WATCH/EURECA instruments were taken into account. The 
failure of unit number 2 on board GRANAT, and the limited 
field of view and the Earth blockage of WATCH/EURECA, 
made it that none of the experiments covered uniformly the 
sky. The WATCH/GRANAT map shows larger exposures to- 
wards the Galactic centre whereas the WATCH/EURECA one 
is under exposured towards the equatorial poles (Brandt 1994, 
Castro-Tirado 1994). If we assume that GRBs occur randomly 
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both in space and time, the probability of detecting a GRB in 
a given direction is proportional to the exposure time spent on 
that region. Therefore for each simulated set of 57 bursts, 45 
of them follow the WATCH/GRANAT exposure map, 10 the 
WATCH/EURECA exposure map and the remaining two bursts 
(representing GRB 920814 and GRB 921022) follow both ex- 
posure maps simultaneously. The simulated WATCH-like sets 
have the same error radii than the real WATCH catalogue. 

2.4. Random and real overlaps 

We call random overlaps to overlaps between the BATSE bursts 
and the simulated WATCH events. The random overlaps pro- 
vide a set of Cij that follows a distribution so-called c ran dom- 
In order to estimate such distribution, the value of the overlap- 
ping functions are calculated for all the overlaps between 
the BATSE sample and 50 simulated WATCH catalogues. It 
shows a mean value < c ran dom > = 0.0098 and a deviation 
Crandom = 0.035. c ran dom provides the expected value of the 
overlapping function when there is a casual overlap between 
two boxes (not due to arise from the same source). The major- 
ity of the random overlaps shows very low values of the over- 
lapping function because they tend to occur at the border of the 
error boxes in the tail of the probability distribution. 

On the other hand, the overlapping function, cy, for each of 
the 27 BATSE- WATCH simultaneous pairs is calculated. The 
distribution of these 27 values of is called e rea i. The mean 
value of the real overlaps, < c rca i >= 0.28932, and the devi- 
ation er rea i = 0.22984. As expected < c rea i > is greater than 
< Crandom >■ This fact can be explained taken into account 
that the probability distributions due to a single GRB detected 
by both experiments tend to be close to each other, compared 
with two GRBs randomly located in the same zone in the sky. 
Thus, the random overlaps tend to occur in the tail of the prob- 
ability distribution, thus forcing cy to be very low. Moreover, 
the lower sensitivity of WATCH in comparison to BATSE im- 
plies that the 27 simultaneous bursts are brighter than the av- 
erage BATSE bursts, (as the radii of the error boxes depend on 
the intensity) and therefore they have smaller error boxes, thus 
making < c rca i > larger than < c ran dom >■ 

The next step is to consider c rca i as the expected distribu- 
tion of for repeaters. This consideration is based on the two 
following assumptions: 

1 . There is little variation with time on the sensitivity of both 
experiments. A change in the sensitivity imply into differ- 
ences in the sizes of error boxes and thus in the Cy values. 
If more accuracy is desirable, then it is necessary to know 
how the sensitivity of both instruments evolves, in order to 
correct the sizes of the error boxes depending on the date 
of detection. 

2. The intensities of different bursts from a repeater source 
do not change significantly in time. Therefore the sizes of 
the repeater error boxes remain approximately the same. A 
more complicated study would deal with the time evolution 
of the repeater sources. 



2.5. Quantification of the number of repeaters 

The set of 1500 Cj's calculated using mock WATCH cata- 
logues follows a Gaussian probability distribution (hereafter 
called So, see Fig. ||). The simulated WATCH catalogues were 
generated only using the exposure maps and they only contain 
accidental overlaps, because the simultaneous GRBs were ex- 
cluded from the sample. Therefore So gives us the expected 
value of the "added correlation" when WATCH and BATSE 
catalogues do not share any source. Assuming that c re al rep- 
resents the expected value of the overlapping function for re- 
peater sources, we can introduce trial repeaters and construct 
the Cj's probability distributions Sn for different number of 
repeaters, N, by the following symbolic expression : 

N 

Sn = So + S ^ (c rca l — C ran dom) 
i=l 

N being the number of repeaters. If we take any "added 
correlation" Cj of N = repeaters with a probability given 
by Sq, and then we add the contribution to Cj of any repeater 
with overlapping function given by c rea i and subtract the con- 
tribution of any random overlap given by c ran dom, we get a new 
"added correlation" Cj . This process can be repeated by intro- 
ducing other real and random overlaps, providing a new set of 
Cj's. Once a trial repeater has been introduced, this set of Cj's 
will follow a different probability distribution from So, called 
Si. Thus, Si provides the expected values of Cj when BATSE 
and WATCH share one source. Similarly, this method can be 
applied for N — {2, 3, ...9} repeaters, in order to obtain the 
distributions of the "added correlations" for different number 
of repeaters, Sn, N — {2, 3, ...9}. Fig. [| shows the probabil- 
ity distributions Sat, N = {0, 1, 2, 3, ...9} obtained using this 
procedure. 

The intersection between Cbw an d the distributions Sn 
provides 9 values called Pn, N — {0, 1, 2. ..9}. Based on the 9 
values of Pn we can obtain the distribution of Pn for AT > 9. 
Taking into account that the maximum number of allowed co- 
incidences is 57 x 1878, the distribution of Pn's can be nor- 

^ 57 x 1878 

malized by imposing 2~2n=i ^ n = T Then, Pn provides 
the probability that the BATSE and WATCH catalogues share 
N sources (see Fig. Q). 

3. Results and discussion 

As it can be seen in Fig. |[ the mean value of the "added corre- 
lation" for the simulated catalogues, < S >— 16.40 ± 1.47, 
is even larger than the "added correlation" for the real WATCH 
and BATSE catalogues, namely Cbw = 13.72. This implies 
that our results agree qualitatively with the absence of common 
sources. 

The fact that we have preferred to simulate WATCH cat- 
alogues instead of BATSE ones is only due to the computing 
time, because it is more efficient to simulate sets of 57 bursts in 
comparison with groups of 1906 members. In spite of this fact, 
the roles of both catalogues were exchanged in order validate 
the method, applying the process explained in section § 2 to 50 
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Fig. 2. The values of the "added correlation" Cj, for the simu- 
lated catalogues j = {1, 2, ...1500}. The solid line represents 
< S >, the mean value of the "added correlation" for the sim- 
ulated catalogues, the long-dashed lines are the ±lcr limits. As 
it is clearly seen the real value of the "added correlation" Cbw 
(represented by the square) is below the ler limit. Our results 
are not compatible with the presence of common sources, as 
expected from the graph. 
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Fig. 3. The Gaussian-like curves represent the probability dis- 
tributions Sn, N = {0, 1, 2. ..9}, and the vertical dashed line 
shows the value of the BATSE- WATCH "added correlation". 
The intersection of Cbw with the probability distributions Sn 
provides the set P N ,N = {0, 1,2. ..9}. 
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Fig. 4. The values of Pn for different number of repeaters. The 
probability of having N repeaters reaches its maximum at N = 
(no common sources) and decreases rapidly with the number 
of repeaters. 

Table 1. The first and third column represent the number of re- 
peaters, N. The second and the fourth ones give the probability 
of having N repeaters. 



N 


Pn (%) 


N 


Pn (%) 





24.8 ±4.0 


5 


6.0 ±2.0 


1 


18.7 ±3.0 


6 


4.5 ±1.7 


2 


14.0 ±2.7 


7 


3.4 ±1.5 


3 


10.5 ±2.5 


8 


2.5 ±1.1 


4 


7.9 ±2.2 


9 


1.9 ±0.9 



not share any source. Thus, our results support the lack of com- 
mon sources. Furthermore, the number of common sources is 
< 9 with a 94% confidence level (see Table 1), which means 
a 15.8% of the whole sample. This percentage is similar to the 
20% upper limit imposed to the IB catalogue (Strohmayer et al. 
1994). The results are also in good agreement with the studies 
carried out with the BATSE 3B (Tegmark et al. 1996) and 4B 
catalogues (Hakkila et al. 1998), which did not find evidence 
of repetition. A possible reason to explain our results could be 
due to the different sensitivity of the experiments, as WATCH 
is sampling the strongest bursts and BATSE is also detecting a 
fainter population. The different populations of objects found 
inside WATCH and BATSE error boxes could support this idea 
(Gorosabel and Castro-Tirado 1998a, 1998b). 



BATSE simulated catalogues. Only with 50 catalogues the val- 
ues obtained for Cbw and < S > differ by less than 5% from 
those obtained when WATCH catalogues were simulated. 

The probability distributions Sn are shown in Fig. [}] and 
the deduced values of Pn are given in Table 1 and displayed 
in Fig. As it is shown in Table 1, Pn decreases with N, 
showing the maximum value when BATSE and WATCH do 



4. Conclusion 

In this study we have developed a method that allows us to 
search for GRBs common to two catalogues of sources, each 
one based on a different instrument. The method makes use 
of the GRBs detected simultaneously by both experiments, so 
it is necessary that the experiments overlap in time. We have 
applied the method to the WATCH (WATCH/GRANAT + 



6 



Gorosabel et al.: Search for common WATCH-BATSE GRBs 



WATCH/EURECA) and BATSE (BATSE 4B + bursts detected 
after August 1996) catalogues. 

We conclude that there is no evidence of recurrent activity 
of WATCH bursts in the BATSE sample. We claim (with a 94% 
confidence level) that no more than a 15.8% of the 57 GRBs 
detected by WATCH are present in the sample of 1905 BATSE 
bursts (excluding the simultaneous bursts). However, the pos- 
sibility of finding repeaters in each single catalogue cannot be 
ruled out. Our results support models which do not predict rep- 
etitions of GRBs (for instance the merging of neutron stars at 
cosmological distances). 
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